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Abstract. Theories of planet formation predict the birth of giant planets in the inner, dense, and 
gas-rich regions of the circumstellar disks around young stars. These are the regions from which 
strong CO emission is expected. Observations have so far been unable to confirm the presence 
of planets caught in formation. We have developed a novel method to detect a giant planet still 
embedded in a circumstellar disk by the distortions of the CO molecular line profiles emerging 
from the protoplanetary disk's surface. The method is based on the fact that a giant planet 
significantly perturbs the gas velocity flow in addition to distorting the disk surface density. 
We have calculated the emerging molecular line profiles by combining hydrodynamical models 
with semianalytic radiative transfer calculations. Our results have shown that a giant Jupiter- 
like planet can be detected using contemporary or future high-resolution near-IR spectrographs 
such as VLT/CRIRES or ELT/METIS. We have also studied the effects of binarity on disk 
perturbations. The most interesting results have been found for eccentric circumprimary disks 
in mid-separation binaries, for which the disk eccentricity - detectable from the asymmetric line 
profiles - arises from the gravitational effects of the companion star. Our detailed simulations 
shed new light on how to constrain the disk kinematical state as well as its eccentricity profile. 
Recent findings by independent groups have shown that core-accretion is severely affected by 
disk eccentricity, hence detection of an eccentric protoplanetary disk in a young binary system 
would further constrain planet formation theories. 
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1. Introduction 

It is well known that circular Keplerian protoplanetary disks are expected to produce 
symmetric double-peaked molecular line profiles (Home & Marsh 1986). Contrary to this 
simple symmetric disk assumption, asymmetric CO line profiles in the fundamental band 
have been observed in several cases. Grid-based numerical simulations of Kley & Dirksen 
(2006) have shown that local disk eccentricity might form in planet bearing disk near the 
gap. The theory of resonant excitation mechanisms in accretion disks of Lubow (1991) 
predicts that the circumstellar disks of close-separation young binaries might become 
fully eccentric due to the orbiting companion. Home (1995) has shown that supersonic 
turbulence might cause observable line profile distortions. As the orbital velocity of gas 
parcels is highly supersonic in accretion disks, the disk eccentricity results in supersonic 
velocity distortions. Therefore distorted molecular line profiles are expected to form in 
giant planet bearing locally eccentric protoplanetary disks (Regaly et al. 2010) and fully 
eccentric circumstellar disks of close-separation young binaries. 
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2. Spectral calculations combined with hydrodynamic simulations 

In order to calculate the CO spectra, we need to know the temperature distribution in the 
disk. In our study the double layer flaring disk model of Chiang & Goldreich (1997) is as- 
sumed. In this model the disk is heated by the stellar irradiation and accretion processes. 
The incident stellar flux heats the disk atmosphere, which reprocesses the stellar light 
and irradiates the disk interior. The accretion processes heats the disk interior directly. 
In this way an optically thick interior and an optically thin atmosphere develop (Fig. 1). 
For low accretion rate, temperature inversion forms, resulting in emission spectra. 
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Figure 1. The double layer disk model in which the disk atmosphere is superheated. 



The fundamental band CO emission is formed in the optically thin disk atmosphere 
above the optically thick cooler disk interior. In the optically thin approximation the 
/(z/, <j)) monochromatic intensity of the radiation at frequency v from a given gas parcel 
at i?, (j) cylindrical coordinates is the sum of the dust continuum of the disk interior and 
the optically thin atmospheric CO emission 

I(v, R, <P, i) = B(u, T d {R))e~ T ^ R ^ + B(u, T g (R))(l - e~ T ^ R ^), (2.1) 

where, B{y, T) is the Planck function, Td(R) and T g (R) are the dust and gas temperatures 
in the disk interior and the disk atmosphere, respectively. Details of the spectral model 
can be found in Regaly et al. (2010). For a disk perturbed by a companion (a planet or a 
secondary star) the velocity field is no longer circular Keplerian, in which case the local 
Doppler shift that affects the atmospheric optical depth (r(i/,R,(p,i)) can be given by 

Av(R, = ^ sin(i) {[u R (R, (j)) cos(0) - u^R, (j)) sin(0)] 

+ [u R (R, 4>) sin(0) + u^R, </>) cos(0)]} , (2.2) 

where u R (R,(j)) and u^{R^4>) are the radial and the azimuthal velocity components of 
the orbiting gas parcels. In order to calculate the local Doppler shift in a perturbed disks 
with embedded planet or a secondary star, we solved the continuity and Navier-Stokes 
equations by the 2D grid-based hydrodynamic code FARGO (Masset 2000). In this way, 
the u R (R,(j)) and u ( f ) (R : <p) velocity components of gas parcels, were provided by the 
hydrodynamic simulations. For simplicity we used a locally isothermal equation of state 
for the gas. The a- type disk viscosity was applied (Shakura & Sunyaev 1987). The disk 
self-gravitation was neglected as the Toomre parameter is well above 1. For more details 
on the hydrodynamic simulations see Regaly et al. (2010). 
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3. Planetary signal: local disk eccentricity due to a giant planet 

According to our simulations the density distribution shows permanent elliptic shape 
near the gap (Fig. [5J left panel), confirming the results of Kley & Dirksen (2006). We 
found that the orbits of gas parcels are eccentric in the vicinity of the planet (Fig. [5J 
middle panel). As the magnitude of the velocity distortion of the gas parcels exceeds the 
local sound speed, distorted non-symmetric line profiles are expected to form. 

We have calculated the fundamental band P10 non-blended CO line for 20, 40 and 60 
degree of disk inclination in a perturbed disk assuming an 8Mj up mass planet orbiting 
1 Mq mass star at 1 AU (Fig. [2j right panel). The line shape is distorted and the relative 
strength of the distortions are increasing with increasing inclination angle. Calculating 
the line profiles with different orbital positions of the planet, we have found that the line 
shapes vary as the planet is orbiting. 

In order to investigate the dependence of the line profile distortions on the model 
parameters, we have calculated the CO emission in several models assuming different 
planetary and stellar masses, and orbital distance of the planet. We found that the 
planetary and the stellar mass affect the line profiles by the same means: the larger the 
planetary or stellar mass, the stronger the distortion. The line profile distortion is found 
to be weakening with increasing orbital distance. According to our calculations, the signal 
of a Jovian planet orbiting a Solar mass star at a distance of 1 AU can be detected with 
CRIRES (Regaly et al., 2010). 




Figure 2. Surface density (left panel) and radial component of velocity (middle panel) distri- 
bution in a planet bearing protoplanetary disk. Distorted V = 1 — > 0P10 CO lines (right panel) 
emerging from the gravitationally distorted disk, assuming 20°, 40° and 60° disk inclinations. 



4. Asymmetric lines: fully eccentric disk due to a secondary star 

Circumprimary disks in young binaries become fully eccentric via the gravitational per- 
turbations of the secondary, in form of resonant excitation mechanisms (Lubow 1991, 
Kley et al. 2008). Thebault et al. (2006) and Paardekooper et al. (2008) have shown that 
the planetesimal accretion might be inhibited by the disk eccentricity in core-accretion 
models. Recently, Zsom et al. (2010) have shown that the dust coagulation process is 
also affected by the disk eccentricity in the core-accretion scenario. 

According to our simulations, the average disk eccentricity is ~ 0.3 independently on 
the binary mass ratio and the magnitude of viscosity, assuming h = 0.05 aspect ratio for 
the disk, which is a reasonable assumption for a circumstellar disk around an average 
TTauri star (Fig. [3j left panel). However, smaller disk eccentricity was found for thick 
(h ^ 0.75) and thin (h ^ 0.03) disks. We have found that the orbital eccentricity of the 
binary above 0.2 inhibits the development of disk eccentricity. 

Calculating the V = 1 — »• 0P10 non-blended CO lines with 20° and 60° disk inclination 
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angle, we have found a maximum ~ 20% red-blue peak asymmetry (Fig.[3j middle panel). 
The profile asymmetry is due to the fact that the inner disk, where the fundamental band 
CO emission is formed, is fully eccentric. The peak asymmetry is inverted twice during 
a full disk precession, which takes about ~ 6.5Pbi n - Moreover, as the disk eccentricity 
slightly varies during one binary orbit, line profile variations might be observed in the 
line wings with ~ Pbin- 
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Figure 3. The circumprimary disk becomes eccentric due to the perturbations from the sec- 
ondary. Asymmetric V = 1 —> 0P10 CO line (middle panel), assuming 20° and 60° disk inclina- 
tions. Trailed spectra (right panel) covering a full disk precession shows the periodic variations 
in asymmetry with ~ 6.5Pbin and line wings with ~ Pbin- 



5. Summary 

High-resolution 4.7/im CO molecular spectra provide us a tool to reveal disk eccentricity 
within the 3 AU regions, where terrestrial planets are expected to form. We have shown 
that Jovian planets still embedded in their protoplanetary disk can indirectly be detected 
with the 4.7/im CO line profile distortions caused by the local disk eccentricity. We have 
also shown that significant CO line profile asymmetry and variations are expected from 
fully eccentric disks of close separation young binaries. By revealing disk eccentricity in 
close-separation young binaries the core-accretion model might be constrained further, 
since the planet formation seems to be strongly affected by the disk eccentricity. 
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